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ABSTRACT 

We show how the semiclassical Langevin method can be extended to calculations of higher-than-second cumulants 
of noise. These cumulants are affected by indirect correlations between the fluctuations, which may be considered 
as " noise of noise." We formulate simple diagrammatic rules for calculating the higher cumulants and apply them 
to mesoscopic diffusive contacts and chaotic cavities. As one of the application of the method, we analyze the 
frequency dependence of the third cumulant of current in these systems and show that it contains additional 
peculiarities as compared to the second cumulant. The effects of environmental feedback in measurements of the 
third cumulant arc also discussed in terms of this method. 
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1. INTRODUCTION 

Since the nonequilibrium noise in mesoscopic systems became a hot topic in late 80s, the semiclassical Boltzmann- 
Langevin approach has been successfully competing with quantum- mechanical methods in these studies. 1 While 
restricted to certain class of systems and parameters, this approach is generally more simple and physically 
appealing. It can easily handle some of the problems that are extremely difficult to solve by directly using 
quantum-mechanical equations. 

Higher-order correlations of current became a subject of interest for theorists since early nineties. This work 
was pioneered by Levitov and Lesovik, 2 who discovered that the charge transmitted through a single-channel 
quantum contact obeys a binomial distribution. Based on their quantum-mechanical formulas, higher cumulants 
of current were calculated for a variety of multichannel phase-coherent systems. 3-5 Meanwhile a wide class of 
these systems, like diffusive conductors and chaotic cavities, allows a semiclassical description of their average 
transport properties and second cumulant of noise. 6, 7 Hence it is of interest to have a fully semiclassical theory 
for higher cumulants of noise in these systems. 

Below we present an extension of the Boltzmann-Langevin approach that allows one to calculate the statistical 
properties of noise beyond the second cumulant. This extension is based (a) on the time-scale separation between 
the correlation time of microscopic sources of noise acting upon the system (" noise generators" ) and the response 
of the system to them and (b) on the smallness of fluctuations. It has been initially proposed for the higher 
cumulants of current in mesoscopic diffusive wires, 8 where the resulting terms could be directly mapped on the 
quantum-mechanical diagrams for these quantities. 9 Shortly after this, the cascade expansion was postulated for 
a chaotic cavity with quantum contacts and tested against quantum-mechanical results for the third and fourth 
cumulants of current. 10 Very recently, a general proof of this expansion was given for semiclassical cumulants of 
arbitrary order. 11 ' 12 

An application of the semiclassical method to calculations of the frequency dependence of the third cumulant 
of current noise in mesoscopic diffusive wires 13 and chaotic cavities 14 gave very unexpected results. Unlike the 
conventional noise and ac electric response that exhibit a dispersion only at the inverse RC time of the system, 
the third cumulant of current has also a dispersion at the inverse dwell time of an electron on the system. 
This dispersion is due to slow fluctuations of the distribution function that do not violate electroneutrality and 
are akin to fluctuations of local temperature. These fluctuations do not directly contribute to the current and 
therefore are not seen in conventional noise, but they modulate the intensity of noise sources and therefore 
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manifest themselves in higher correlations of current. As the dwell time is generally much larger than the RC 
time, the resulting dispersion may be experimentally observed. 

The cascaded expansion of higher cumulants is valid not only for systems described by Boltzmann equation 
but also for other systems that satisfy the conditions of time-scale separation and smallness of fluctuations. 
Recently, this method has been used for calculating the feedback of the external circuit in measurements of the 
third cumulant of an arbitrary mesoscopic resistor. 15 

The paper is organized as follows. In Section 2 we describe the general cascaded formalism. We give a brief 
overview of the standard Boltzmann -Langevin method, then we derive the expression for the third cumulant using 
a generalized Fokker-Planck equation. After this, we formulate the general diagrammatic rules for calculating 
cumulants of arbitrary order. In Section 3 we calculate the frequency dependence of the third cumulant for 
a mesoscopic diffusive wire and a chaotic cavity. Finally, we address the effect of the external circuit on the 
measurements of the third cumulant. The main conclusions are summarized in Section 4. 



2. GENERAL FORMALISM 

The basis for the semiclassical description of kinetics is the existence of two well separated time scales, one of 
which describes a "slow" classical evolution of the system and the other describes "fast" quantum processes. For 
example, the collision integral in the Boltzmann equation may be written as local in time because quantum- 
mechanical scattering is assumed to be fast as compared to the evolution of the distribution function. For this 
reason, it is possible to describe the noise associated with electron scattering by Langevin sources that are 5- 
correlated in time. On the other hand, using the semiclassical kinetic equation implies coarse-grain averaging 
over N = (ApAx/2irh) 3 states located in an elementary cell of phase space whose characteristic dimensions Ap 
and Ax are small as compared to the characteristic scales of the problem. This suggests that the cumulants of 
the distribution function / should decay with number as l/N™^ 1 . 

2.1. Boltzmann Langevin Method 

The second cumulants of noise are conveniently described by the Boltzmann-Langcvin equation, 6 which is ob- 
tained by linearizing the standard Boltzmann equation with respect to a fluctuation of the distribution function 
5f. It reads 

6f(p,r,t) + 5I=-eSEv^ + J ext , (1) 

where SI is the linearized collision integral and J ext is the Langevin source that accounts for the randomness of 
electron scattering. This equation should be supplemented by a self-consistency equation for the electric field 

V6V = 4tt8 P , 5p(r,t) = eJ2Sf(p,r,t). (2) 

p 

The key issue in the Boltzmann-Langcvin method is a derivation of the correlation function of Langevin sources. 
Several decades ago,Kogan and Shulman 6 proposed a physically appealing expression fort this function. Because 
of a short duration of scattering events, the Langevin sources are (5-correlated in time. Furthermore, since 
different scattering events are independent, the scattering between each pair of states (p,r) and (p',r) presents 
a Poissonian process whose cumulants of any order are proportional to its average rate 

J(p - p') = VF(p,p')/(p,r,t)[l - /(p',r,t)]. (3) 

Hence the cumulants of extraneous sources of the corresponding order related with randomness of scattering may 
be written as the sums of incoming and outgoing scattering fluxes taken with appropriate signs. For example, 
the second cumulant is given by 6 

(( J ext ( Pl , n, h) J ext (p 2 , r 2 , t 2 ))) = 5( ri - r 2 )S(h - t 2 ) | S P1P2 £[J(pi - p') + J(p' - pi)] 

I p' 



d_ + v _9 +eEv _3 

dt dr de 



-J(pi-»p 2 )-J(P2-»Pi)|. (4) 

To calculate, e.g., the correlator of current fluctuations (5j(ri, ti)Sj(r 2 ,t 2 )), one has to take two formal solutions 
of Eq. (1) in 5f, calculate the two fluctuations of current via 

5j(r,t)=e^v5f(p,r,t), (5) 
p 

multiply them, and average the product using Eq.(4). As the fluxes (3) depend on the average distribution 
function,the resulting noise is also a functional of /. 

The Boltzmann-Langevin method has been very successful in predicting the 1/3 suppression of shot noise in 
mesoscopic diffusive contacts 16 and analyzing the effects of electron-electron scattering on it. 17 

2.2. Calculation of the Third Cumulant 

One might think that all higher cumulants of order n could be treated similarly just by averaging the product 
of n formal solutions of Eq. (1) with the corresponding analog of Eq. (4). However this is not the case. The 
problem is that the correlator of Langevin sources (4) as well as higher-order correlators of 5 J ext are themselves 
functions of /, which fluctuates in time, and this results in additional correlations between <5/'s. To elucidate 
this point, we consider a simple example of a Markov process where a random quantity x(t) is described by a 
Langevin equation 

x = --(x-x) +£(t) (6) 

T 

and £(t) is a random extraneous force with zero mean (£(t)) = 0. Suppose that an nth cumulant of £ is inversely 
proportional to a large number TV™ -1 and explicitly depends on x: 

<£(ti)£(t 2 )> - ^5(h - t 2 )Q( Xl ), (£(ti)£(t 2 )£(t 3 )> - ^6(h - t 2 )5(t 2 - t 3 )R( Xl ). (7) 

This suggests that as At — > 0, the cumulants of increment of x are proportional to At and equal 

(x(t+At)-x(t)) = -^[ X (t)-x], ([x(t+At)-x(t)} 2 ) = ^Q(x(t)), ([x(t + At)-x(t)} 3 ) = ^R(x(t)). 

It's a textbook knowledge that the probability w(x,t) in this case obeys a generalized Fokker-Planck equation 
in the form 

dw _ 1 d 1 d 2 {Qw) 1 d 3 (Rw) 

lx~Tlte [{x ~ x > wi + w~dx~ 2 6iv 2 dx 3 ■ ( ) 

By multiplying both parts of Eq. 8 by 5x = x(t) — x, Sx 2 , and Sx 3 and integrating them over x by parts the 
corresponding number of times, one obtains equations of motion for the conditional cumulants of x 

Y x + • • •) - o, (± + 1) (Sx 2 (t) ...) = l(g W ), (9 ) 

and 

|_ + ^ {5x 3 (i) ...) = ^(Sx(t)Q(t)) + ±(R(t)), (10) 

where . . . stand for any functions of x taken at previous instants of time, (Q(t)) = (Q(x + Sx(t))), and (R(t)} = 
(R(x + Sx(t))). By subsequently solving the equations for lower cumulants with senior cumulants as the initial 
conditions, one obtains the expressions for the two- and three-time correlation functions of Sx for t\ > t 2 > t 3 
in the form 

t 

{SxitJSxfa)) = ^ J dt' K{t 1 -t')K{t 2 -t'){Q{t')) (11) 



and 



ti 

(5x{t x )6x{ti)5x<t z )) = ^ J dtfKiti-^Kfo-WQMSxfa)) 

t 3 

+ J dt' K(t 1 -t')K(t 2 -t')K(t 3 -t') 



l(Q(t')6x(t')) + ^{R(t')) 



where K(t - t') = exp[-(i - t')/r]. 

Expand now Q(t') and R(t') in powers of Sx(t') about x. By substituting these expansions 



(12) 



Q(t) = Q(x) + 5x(t') ^ 



+ . . . , R(t) = R(x) + 5x(t') ^ 



+ 



into Eqs. (11) and (12), one obtains a system of equations where cumulants of all orders are coupled together 
and which cannot be solved in a general case. However one can make use of the smallness of l/N. Retaining 
only leading terms of the order of l/N n ~ 1 in nth cumulants, it is possible to solve the equations recursively 
by expressing the higher cumulants in terms of the lower ones. Thus the expression for the second cumulant 
becomes just 



(6x(h)5x(t 2 )) 



1 J dt 1 K(h-l/)K(t 2 -lf)Q(x). 



(13) 



This is exactly what the standard Langcvin approach with Q = Q(x) gives. 

The expression for the third cumulant appears to be more involved. Taking into account that K (ti 
t')(Q(t')5x(t')) = (5x(ti)Q(t')) if U > t', it is easily brought to a form 



12 

(Sx^Sx^Sx^s)) = jj J dt' K(h - t')K(t 2 - t') 



„ dQ 
dx 



(5x(t')Sx(t 3 )) 



t 3 t 2 
^ J dt' K(t 2 - t')K(t 3 - t') ^ (fe(i')fe(ii)) + ^ j dt' K(h - t')K(t 3 - t') ^ 



(6x(t')Sx(t 2 )) 



13 

+ J dt' K(h-lf)K(t2-1?)K(t3-1?)R(x). (14) 

— oo 

We note that the first three terms of Eq. 14 present convolutions of the functional derivatives of the second 
cumulant with another second cumulant, so that symbolically, this equation may be recast in a form 

(SxxSx^x^ = (SxxSx^x^q^ + P 123 | S ^ X ^ X2 ^ (Jx<5x 3 ) j , (15) 

where P\ 23 denotes a summation over all inequivalent permutations of indices (123) and S(...)/Sx denotes 
a functional derivative with respect to x. The first term in Eq. (15) presents a direct contribution to the 
third cumulant of x from the third cumulant of extraneous sources, and the second term presents the cascade 
corrections, which result from modulation of the second cumulant by fluctuations of x. 

If x denotes a semiclassical distribution function /(p, r, t), one may roughly imagine these cascade corrections 
as fluctuations of Nyquist noise of a resistor caused by fluctuations of its temperature. It should be noted that 
the conditions for validity of the cascaded approach are the same as the conditions of validity of Boltzmann 
description itself and no additional assumptions are needed. 
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Figure 1. The second cumulant and the two contributions to the third cumulant of current. The external ends correspond 
to current fluctuations at different moments of time and the dashed lines, to fluctuations of the distribution function. The 
full circle and triangle correspond to cumulants of extraneous currents and the empty triangle to the functional derivative 
of the second cumulant. 

2.3. Diagrammatic Expansions for Higher Cumulants 

Similar expressions are valid for higher cumulants. Apart from direct products of n formal solutions of the 
Langevin equation, nth cumulant also contains cascade corrections from an interplay between lower cumulants. 
The cascade corrections are conveniently presented in a diagrammatic form (see Figs. 1 and 2). All diagrams 
present graphs whose outer vertices correspond to different instances of x and whose inner vertices correspond 
either to cumulants of extraneous currents or their functional derivatives. The number of arrows outgoing from 
an inner vertex corresponds to the order of the cumulant and the number of incoming arrows corresponds to the 
order of a functional derivative. Apparently, the difference between the total order of cumulants involved and 
the total number of functional differentiations should be equal to the order of the cumulant being calculated. 
As there should be no backaction of higher cumulants on lower cumulants, all diagrams are singly connected. 
Therefore any diagram for the nth cumulant of the current may be obtained from a diagram of order m < n by 
combining it with a diagram of order n — m + 1, i.e., by inserting one of its outer vertices into one of the inner 
vertices of the latter. Hence the most convenient way to draw diagrams for a cumulant of a given order is to 
start with diagrams of lower order and to consider all their inequivalent combinations that give diagrams of the 
desired order. The analytical expressions corresponding to each diagram contain numerical prefactors equal to 
the numbers of inequivalent permutations of the outer vertices. 

At first, these rules had been postulated, 8 ' 10 but very recently they were proved independently by Jordan et 
al. using the stochastic path integral approach and by Gutman et al. 12 using the supersymmetry method. 



3. APPLICATIONS OF THE METHOD 

We now concentrate on the two most important mesoscopic systems that allow a semiclassical treatment, i.e. 
diffusive wires and classical chaotic cavities. Both these types of multichannel systems have a complicated 
internal dynamics characterized by essentially different time scales. In the case of a diffusive wire, thire is a large 
separation between the elastic scattering time and the time of diffusion of an electron across the wire. In the 
case of a chaotic cavity, there is a large separation between the time of flight through the cavity and the dwell 
time of the electron on it. The fluctuations in these systems are relatively small because of a large number of 
quantum channels in them, so the cascaded Langevin approach may be conveniently applied to them. 

Of a primary interest to us will be the frequency dependence of the third cumulant of current 

P 3 {U 1 ,U 2 ) = J d(h-t 2 ) J d{t 2 -t 3 ) explicit! -t3) + iU2{t2-t3)](SI{t 1 )5I(t2)SI{t 3 )) (16) 

in these systems. The point is that in the case of a good conductor with a short screening length, the dwell 
time of an electron in the system is seen neither in its linear response nor in the second cumulant of current 
because the charge transport is primarily controlled by the screening effects whose characteristic time scale is 
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Figure 2. The contributions to the fourth cumulant of the current. Dashed lines correspond to fluctuations of the distri- 
bution function. Full circles, triangles and squares correspond to the second, third, and fourth cumulants of extraneous 
currents. The empty triangles and squares present their functional derivatives. 

the RC time of the system. Typically, this time is so short that the resulting dispersion of noise cannot be 
experimentally observed. In contrast to this, the third cumulant of current is affected not only by fluctuations 
of the total electron number at a given point, but also by fluctuations of their distribution in energy, which may 
affect the intensity of noise sources. The latter fluctuations have a much longer relaxation time since they do not 
violate the electroneutrality and therefore may result in a dispersion at experimentally accessible frequencies. 

3.1. Diffusive wires 

Consider a quasi-one-dimensional diffusive wire of length L and conductivity a with a constant voltage drop V 
across it. All dimensions are assumed to be much larger than the elastic mean free path and the screening length 
in the metal. We also assume that dimensions are much smaller than the electron-phonon scattering length, but 
take into account electron-electron scattering. It is assumed that the contact is located close to an external gate, 
so it is described by a single capacitance. In the case of a diffusive metal, the distribution function f(s, r) and 
its fluctuations are almost isotropic in the momentum space, and the Boltzmann-Langevin equation reduces to 
a stochastic diffusion equation 



(J- t - DV 2 ^j Sf - SI ee = -eS4>^ - V<5F""P - SF ee , 



(17) 



where D is the diffusion coefficient, SI ee is the linearized electron-electron collision integral, and SF lmp and 
SF ee are random extraneous sources associated with electron-impurity and electron-electron scattering. The 
fluctuation of the electric potential 5(j) that appears in this equation should be calculated self-consistently from 
the Poisson equation 

V 2 5(t> = -Air5p, (18) 
where the fluctuation of charge density Sp is given by 



Sp = eN F (^J ds6f(, 



e) + eh 



(19) 



and where Np is the Fermi density of states. In the case of a quasi-one-dimensional contact, a solution of Eqs. 
(17) - (19) is of the form 18 



84>{x,w) = ^ (V 2 + 10JRC/L 2 )- 1 A J rfVjjr f (r), 



(20) 



where x is the coordinate along the contact, So is the cross section area of the contact, C and R are the 
capacitance and the resistance of the contact, and 



Sj ext — eNp J deSF imp . (21) 



A fluctuation of the total current at the left end of the contact equals 



51 = a d 2 r ± 



2 d6< j>(x,w) 
dx 



(22) 

x=-L/2 



Making use of the correlation function of extraneous sources 



(SFr P (e, r)SF; mp (e', r')) w = 2^-6(r r')6(s s')6 af3 f(s, r)[l - f(e, r)], (23) 

one easily obtains the second-order correlation function for the fluctuations of the current as a functional of the 
distribution function f. 16 At frequencies much smaller than (RC)^ 1 it is of the form 

(J/(wi)<5/(w2)> = A-k5(w 1 +w 2 ){RL)- 1 / dx [ dsf{e,x)[l-f(s,x)}. (24) 



To calculate the third cumulant of current in a diffusive wire, one may use the cascaded Langevin approach 
described in the previous section with / in place of x. In the case of diffusive metal, the direct contributions from 
the higher-than-second cumulants of Langevin sources to the corresponding cumulants of current are negligibly 
small because these cumulants are proportional to the inverse elastic scattering time r _1 (see Section 2.1) and 
each solution of the Boltzmann-Langevin equation for the anisotropic part of Sf gives an additional factor of 
r, so that the nth cumulant of transport current would be proportional to t"^ 1 instead of r. Hence all higher 
cumulants of current for diffusive contacts are dominated by diagrams constructed of the second cumulant of 
Langevin sources and its functional derivatives. As has been shown in Section 2.2, the third cumulant in this 
case is of the form 

(6I(t 1 )6I(t2)5I{t 3 ))=P 1 2 3 y dtjdej d\ (6f(e,r,t)6I(t 3 ))} . (25) 

We restrict ourselves to the two limiting cases of absolutely elastic scattering and of a strong electron-electron 
scattering. In the case of purely elastic scattering, SI ee and Sf ee in Eq. (17) may be dropped. Furthermore, 
we consider frequencies much smaller than (i?C) _1 , so the frequency-dependent term in Eq. (20) may be 
neglected and the solution of these equations is straightforward. Using the well-known expression for the average 
distribution function 

= (\ + t) Me + eV/2) + (\-^\ fo(e - eV/2), (26) 



2 LJ J ' ' ' \2 L 
where /o is the equilibrium Fermi distribution, one obtains that in the frequency range of interest, 

^{u 1 ,u 2 )=P(u 1 )+P(u2)+P(-u 1 -u 2 ), (27) 

where 

4 e 2 V gMfiL 2 + 30) sinh(g^) - 8(^£ 2 + 6) cosh^L) + 2^£ 2 + 48 
P rf (o,)--- — ———— (28) 

at eV > T and 

4 e 2 V 2cosh(q LU L) - q u Lsmh(q w L) - 2 

rel\W) — -— „ „ . —. — . 

3 R q^L 6 sinh(g w L) 

at T > eV. Here q u = {iuj/D) 1 ' 2 . At u = these expressions give -(l/45)e 2 V/R and -(l/9)e 2 V/R, which 
corresponds to ^(0, 0) = — (l/15)e 2 V/R and ^(0, 0) = — (l/3)e 2 V/R. These zero-frequency results arc in 
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Figure 3. The real and imaginary parts of the ratio P(u>)/(e 2 I) versus normalized frequency lotd (td = L 2 /D) for 
purely elastic scattering at eV 3> T (solid lines), hot-electron regime at eV S> T (dashed lines), purely elastic scattering at 
eV <C T (dotted lines), and hot-electron regime at eV <C T (dash-dotted lines). Inset: charged and uncharged fluctuations 
of the distribution function /. Charged fluctuations have a short relaxation time and contribute to fluctuations of current 
SI. Uncharged fluctuations do not contribute to SI directly but affect the intensity of noise sources. They decay only via 
slow diffusion and result in the low- frequency dispersion of the third cumulant. 

agreement with the results of Gutman and Gefen. 9 At finite frequency, Eqs. (28) and (29) become complex- 
valued and tend to zero as i/co at lo — > oo. 

In the opposite limit of a strong electron-electron interaction, the distribution function may be assumed to 
have a Fermi shape with a coordinate-dependent temperature T e (x) and electric potential 4>{x) 



1 + exp 



e — e<j>{x) 
T e {x) 



and a fluctuation 5f can be expressed in terms of fluctuations of these quantities 



(30) 



(31) 



A substitution of Eq. (31) into Eq. (25) and integration over the energy readily gives 



L/2 



(32) 



-L/2 



To calculate the correlator in Eq. (32), we have to obtain a Langevin-type equation for 5T e . To this end, we 
multiply Eq. (38) by e and integrate it over e, like it was done when deriving the equation of heat balance. 17 ' 19 
This gives 



(e 2 <W) =-j 



deeV5F v 



(33) 



By solving this equation together with Eq. (20) and making use of the mean effective temperature 



17 



f e (x) = 



IT* 



1/2 



(34) 
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Figure 4. A chaotic cavity capacitively coupled to a gate 



one obtains 



Phot{u) = - 



l2e 2 V 1 



1 - 




) 



(35) 



at high temperatures T 3> eV and 



Phot(uj) = 



12 e 2 V J (wk + nffllMnk + tt/2) - (-l) fc ] 
7T i? ^ (2fc+l)[7r 2 (2/c+l) 2 +iwL 2 /L»] 



(36) 



at low temperatures 3> T, where J and J\ are Bessel functions of order and 1. The corresponding limiting 
values P 3 (0, 0) = -(3/7r 2 )e 2 y/i? at < T and P 3 (0, 0) = -(8/tt 2 - 9/16)e 2 ^/i? at eV > T coincide with the 
results of Gutman and Gefen. 9 

Figure 3 shows that both real and imaginary parts of P have the most pronounced dispersion in the case of a 
high temperature or for strong electron-electron scattering, i.e. when the local distribution function has a nearly 
Fermian shape. This unexpected result is in a sharp contrast with the dispersion of quantum noise, 21 which 
results from sharp singularities in the energy dependence of the distribution function. This difference comes 
from the different symmetry of the second-order correlation functions appearing in Eq. (25) in the elastic case 
at low temperatures and in the hot-electron regime. In the former case, the correlation function (Sf(e,x)SI) LU is 
an odd function of the coordinate, but in the latter case {ST e (x)SI) u} is an even function of x. 

3.2. Chaotic Cavity 

The cascaded Langevin approach may be successfully applied not only to systems described by Langevin equations 
with Poissonian statistics of extraneous sources, which takes place in the case of a Boltzmann kinetics. It can 
be equally well applied to systems with arbitrary statistics of noise sources provided that the conditions of 
time-scale separation and smallness of the fluctuations are observed. An example of a semiclassical system with 
non-Poissonian sources of noise is a classical chaotic cavity, i.e. a metallic island of irregular shape connected to 
the electrodes L, R via two quantum point contacts of conductances Gx ; ij » e 2 /h and arbitrary transparencies 
^l,r (see Fig. 4). As the dwell time of an electron in the cavity t d — c 2 N f /(Gl + Gr) is much larger than the 
time of flight through the cavity, the electrons in the cavity lose memory of their initial phase and are described 
by an energy-dependent distribution function f(e,t). Hence the low- frequency dynamics of / is described by 
a semiclassical equation, while quantum contact serve as generators of noise with essentially non-Poissonian 
statistics. The fluctuations of the electric current in the left and right contacts are given by equations 




(37) 



where (Il)s and (Ir) s are the energy-resolved random extraneous currents generated by the left and right 
contacts. The fluctuation of the distribution function 5f(e) obeys a kinetic equation 




(38) 
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Figure 5. A 3D plot of ReP3(a>i,a>2) for Gl/Gr = 1/2, Yl = Tr = 3/4, tq = 1/3, and td = 10 (dimensionless units). 

where to — e 2 NF/{GL + Gn) is the dwell time of an electron in the cavity, Np is the density of states in it, and 5U 
is a fluctuation of the electric potential of the cavity. This fluctuation is obtained from the charge-conservation 
law 

dSU 1 dSQ 1 f , ,,~ , ,~ . . 

where C is the electrostatic capacitance of the cavity. Equations (37) - (39) suggest that different parts of 5f 
are described by different time scales. The relaxation of electrically neutral fluctuations is described by the 
characteristic time to , whereas the fluctuations of charge are described by tq = [{Gl + Gr) jC + 1/ To] _1 , which 
is much shorter than to for good conductors. 

If the distribution function of electrons in the cavity /(e, t) were not allowed to fluctuate, the contacts would 
be independent generators of current noise whose zero-frequency energy-resolved cumulants ((/£ R )) £ could be 
obtained from a quantum-mechanical formula 

((I2, R ))e = ^ ln{l + V L . R f L . R (e)[l f(e)](e-** 1) 

+T L , R f(e)[l - /L,«(c)](e ex - 1)}| X=0 . (40) 

If the voltage is high enough, the noise of isolated contacts can be considered as white at frequencies at which 
the distribution function / fluctuates. This allows us to consider the contacts as independent generators of white 
noise, whose intensity is determined by the instantaneous distribution function of electrons in the cavity. Based 
on this time-scale separation, we perform a recursive expansion of higher cumulants of current in terms of its 
lower cumulants as described in Section 2. In the low-frequency limit, the expressions for the third and fourth 
cumulants coincide with those obtained by quantum-mechanical methods for arbitrary ratio of conductances 
Gl/Gr and transparencies r^^. 10 Very recently, similar equations for the zero-frequency limit were obtained 
as a saddle-point expansion of a stochastic path integral. 20 Consider now the frequency dependence of the 
third cumulant. We will be interested in the case of a good conductor where the charge-relaxation time tq is 
much shorter than the dwell time to. Unlike the second cumulant of current, the third cumulant /^(wi, W2) in 
general exhibits a strong dispersion at u)\ t 2 ~ l/r^. 14 In this limit, 
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Figure 6. Ps(uj,0) as a function of u> for Gl/Gr = 1, Tl = 1, = 0, and td = 10 (dimensionless units). The solid, 
dashed, and dash-dot curves correspond to tq — 1/2, tq = 3, and to weak electrostatic coupling tq = td (C = oo). 
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This dispersion vanishes for symmetric cavities and cavities with two tunnel or two ballistic contacts. For a 
cavity with two tunnel contacts, the white-noise sources (40) arc linear functionals of the distribution function 
and hence are not affected by charge-neutral fluctuations. For the case of two ballistic contacts, Eq. (40) depends 
only on / and not on /j. Then the low- frequency dispersion does not show up either because fluctuations of 
current and distribution function are uncorrclated due to the symmetry of the noise sources. The shape of 
Pz(iO\, uj 2 ) essentially depends on the parameters of the contacts. In particular, for a cavity with one tunnel and 
one ballistic contact with equal conductances Gl — Gr = G it exhibits a non-monotonic behavior as one goes 
from u>\ — u>2 = to high frequencies. A relatively simple analytical expression for this case may be obtained if 
td S> tq and one of the frequencies is zero: 
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(42) 



The Pz{u, 0) curve shows a clear minimum at to <~ (tdtq) 1//2 and the amplitude of its variation tends to £3(0, 0) 
as t q /t d -> 0. 



3.3. Effects of External Circuit on Measurements of the Third Cumulant 

So far we applied the cascaded method to systems where the noise was governed by a microscopic distribution 
function of electrons. Now we give an example where the noise in a system is determined by a macroscopic 
quantity. 

Recently, Levitov and Rcznikov 22 predicted that the third cumulant of noise in a tunnel contact should present 
a linear function of voltage independent of temperature. However when Reulet et al. 23 tested these predictions 
experimentally, they found very poor agreement of their data with theory. The reason for this description has 
been explained by Beenakker et al. 15 by the feedback of the external circuit using the cascade approach. 

In experiments of Reulet and co-workers, the actual measurable quantity was the voltage drop across a 
macroscopic resistor R s connected in series with a mesoscopic contact with a resistance i? , which produced a 




Figure 7. A mesoscopic contact in series with a macroscopic resistor. The cumulants of current of the contact depend 
on the fluctuations of potential 5<j>. 



non-Gaussian noise (see Fig. 7). Naively, one could expect that the third cumulant of voltage drop would be 
just proportional to the third cumulant of current generated by the contact, but this is not the case. 

Indeed, consider the macroscopic resistor and the contact as independent generators of noise currents Sis and 
SIq. The noise of the macroscopic resistor is Gaussian with the only nonzero second cumulant ((if)) = 2T/R S 
representing its equilibrium noise, and the noise of the contact has nonzero second and third cumulants ((Iq)) 
and ((Iq)), which depend on the voltage across it. The Kirchhoff's law for the fluctuation of the total current in 
the circuit may be written in a form 

W--£+*fc«|* + «b, (43) 
Hs no 

where S(j) is the fluctuation of electric potential to be measured. Hence 

4 *-^<"»-«'»- <44) 

One obtains the second cumulant of voltage just by multiplying and averaging two solutions (44), which gives 

RqRs \ f ,/ T 2\\ , r> T 



Calculations of the third cumulant of voltage according to Eq. (15) give 

W 3 » = - ( «'o 3 » + 3 ( ( 2^ + (d 2 o))) (46) 



Ro + RsJ ' \R + RsJ V Rs 

The first term in this expression is the naively expected result, and the second term presents the cascade 
correction that results from modulations of the second cumulant of current generated by the mesoscopic contact 
by fluctuations of voltage across it. Equation (46) suggests that the feedback of the external circuit is very 
essential even at Rs <C Ro- For a tunnel contact where ((Iq)) = {e<f>/Ro) coth(e^>/2T) and ((Iq)) = e 2 (j)/R , it 
results in a change of slope of (((f> 3 )) versus V curve from negative at eV <C T to positive at eV S> T. This is 
exactly what is observed in experiment. 



4. SUMMARY 

In summary, we demonstrated a semiclassical method for calculating higher cumulants of noise in systems where a 
large separation exists between the characteristic times describing the correlations of microscopic sources of noise 
and the characteristic times describing the dynamics of averages. The calculations of the frequency-dependent 
third cumulants in a number of mesoscopic systems show that their dispersion may be very different from that of 
conventional noise and is very sensitive to fine properties of these systems that almost do not manifest themselves 
otherwise. As the conditions for validity of the method are rather general, it may be applied to macroscopic 
systems, too. 
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